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Abstract 

For the asymmetric simple exclusion process rjt on the integer lattice with 
^ ' two-sided Bernoulli initial condition, we derive exact formulas for the following 

^ . quantities: (1) F{r]t{x) = 1), the probability that site x is occupied at time t; (2) 

^ I the correlation function F{r]t{x) = l,r/o(0) = 1); (3) the distribution function 

'sj- ■ for Qt, the total flux across at time t, and its exponential generating function. 

I. Introduction 



O 
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Since its introduction by Spitzer [TT], the asymmetric simple exclusion process 
^ ' (ASEP) has attracted considerable attention in both the mathematics and physics 

. literature due to the fact it is one of the simplest lattice models describing transport 
far from equilibrium [2l[7l[8]. Recall that ASEP on the integer lattice Z is a continuous 
time Markov process rjt, where (a;) = 1 if a; G Z is occupied at time t, and rit{x) = 
if X is vacant at time t. Particles move on Z according to two rules: (1) A particle at 
X waits an exponential time with parameter one (the particle's "alarm clock"), and 
then chooses y with probability p{x,y)] (2) If y is vacant at that time it moves to 
y, while if y is occupied it remains at x and restarts its clock. The process is called 
"simple" since the jumps are restricted to nearest neighbors: p{x,x + 1) = p and 
p{x,x — 1) = q = 1 — p. All clocks are independent of each other and one need not 
worry that two clocks ring simultaneously since this event has probability zero. A 
rigorous construction of this process can be found in Liggett [7j. We note that the 
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process is called the totally asymmetric simple exclusion process (TASEP) if jumps 
occur only to the left (g = 1) or only to the right {p = 1); and ii p = q = 1/2, the 
process is called the symmetric simple exclusion process (SSEP) |7l|8]. 

In previous work |l3] the authors considered ASEP on the integer lattice Z with 
step Bernoulli initial condition: initially a site in Z+ (the positive integers) is occu- 
pied with probability p while sites in Z~ (the nonpositive integers) are unoccupied. 
Formulas were obtained for ¥{xm{t) < x), the probability that at time t the mth 
particle from the left is at a site < x, with consequent asymptotic results for this 
probability. In this paper we consider the more general situation where initially a site 
in Z+ is occupied with probability p+ while a site in Z~ is occupied with probability 
p_. This initial condition is called the two-sided Bernoulli initial condition. Now it 
makes no sense to speak of the mth particle from the left. Instead we find formulas 
for P(?7j(x) = 1), the probability that site x is occupied at time t, and the correlation 
function F{rit{x) = 1, ?7o(0) = 1). We also derive formulas for the generating function 
(e^*^*), where Qt is the total flux across at time t. 

For the special case of TASEP with two-sided Bernoulli initial condition much is 
already known including various limit theorems [U [10] . (Also see [1] for the stationary 
case p_ = p+.) A simplifying feature of TASEP is that it is a determinantal process 
[6j thus permitting the application of random matrix theory techniques. In |l3] it 
was shown that the limit theorems in [1] for TASEP extend to ASEP for the special 
case of step Bernoulli initial condition (p_ = 0, p+ > 0); and thus, these results [11] 
establish a KPZ universality theorem for ASEP. Whether the TASEP limit theorems 
of [1] for two-sided Bernoulli initial condition can be extended to ASEP from the 
formulas derived in this paper remains to be seen. 

We shall assume throughout that p, q ^ 0. As in [13] we begin with a formula 
derived in [T2| for Py(xm(i) < x), the probablity distribution for the mth particle 
from the left given an initial finite configuration Y; then we average over all F in a 
finite subset of Z; and then we let the subset of Z become unbounded. The formula 
for Py(xm(t) < x) involved multiple integrals over large contours, and this allowed 
Y to become unbounded on the right in [13]. Now, if Y is to become unbounded on 
left and right, we must sum over both large and small contours, so the formula for 
^Y{xmit) < x) must first be recast. This is done in Sec. II. In Sec. Ill we deduce the 
formula for Py(?7j(x) = 1), the probability that site x is occupied at time t given the 
initial configuration Y. In Sec. IV we average over Y and take the limit to obtain 
the formula for ¥{rit{x) = 1) for two-sided Bernoulli initial configuration on Z. The 
formula for the correlation function requires only a small change. Finally, in Sec. V 
we deduce formulas for probability distribution of the total flux across the origin at 
time t and for its exponential generating function. For the special case of SSEP we 
connect this last result with a recent paper of Derrida and Gerschenfeld |3j. 
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II. Finite Configuration, Large and Small Contours 



We begin with Theorem 5.2 of [T2] which gives a formula for Py(xm(t) < x), the 
probability in ASEP when the initial configuration is a (deterministic) finite set Y. 
To state it we introduce some notation, which here will be slightly different. 

First, we set t = p/q and recall that the r-binomial coefficient ['^j^ is defined as 

(1 - r^) (1 - r^-^) • ■ ■ (1 - r^-"+i) 
(1 - r) (1 - r2) . . ■ (1 - r") 

when is a positive integer, 1 when n = (empty products are always defined to 
be 1), and when n is a negative integer. 

We define 



and then 



/(x, fc, = /(x, fc, ei, . . . , a) = n /te, 0) n TTT-^ 



A/(x, k, = /(x, fc, - /(x - 1, fc, = (1 - n C') ^> 0- 

All indices in the products run from 1 to k. Notice that / and A/ depend on t, 
although it is not displayed in the notations. 

Finally, given two sets of integers U and V we define 

(t{U, V) = #{(m, v) : u eU, V eV, and u > v.}. 

(If V = [1, and U C [1, A^] then a{U, V) equals the sum of the elements of U. 
Hence the notation.) 

Theorem 5.2 of [12j states that when g 7^ 0, 

\Y\ 

Py(x„(t) = x) = 5^c„,fc 5^ W^'^) / ■■■/ AI{x,k,Oll^r'd^i---d^k, (1) 



k=l SCY 
|S|=fc 



where 

S = {si, . . . ,Sfc}, 



^This is short for the set {!,..., N}. We use this notation consistently. 
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1 \m k{k~l)/2 m{m—l)/2—km 



k-1 

m — 1 



(2) 



and Cfi is the circle with center and radius R, which is assumed so large that the 
denominators p + qC,iC,j — C,i are nonzero on and outside the contours. The inner sum 
is taken over all subsets S of Y with IS"! = k. Observe that Cm,k = when k < m. 

If we sum (II]) on x from — oo to x, as we may do since i? > 1, we obtain 



Py(x„(t) < x) 



5^ r'^^^'^) / ■■■/ nx,k,Oll^rd^i---d^k. (3) 



k>l SCY 

~ |S|=fe 



Clearly ([3]) implies ([T]) immediately . 

Because of the appearance of the factors J| • in the integrands and because 
R is large we expect to be able to allow Y to be unbounded on the right but not on 
the left. So we want integrals over both large and small contours. Lemma 5.1 of [12] 
gives a formula which expresses an integral over a large contour as a sum of integrals 
over small contours, and we shall be able to apply it here. 

We shall assume throughout that F = y_ UF^ where Y_ is to the left of Y^. Then 
each S* C F is a union S- U 5*+ with 5*^ C Y±. Set kj- = 15*^1, so A; = A;_ + k^. Instead 
of using indices i = 1, . . . , k for each S with 15*1 = k, we shall use /c_ negative indices 
— . . . , — 1 and k+ positive indices 1, . . . , /c+, with Si E S- if i < and Si G 5+ if 
i > 0. In this notation ([3]) becomes 



fyiXmit) <X)= ^ Cm,k ^ T 



a{S-,Y-)+a{S+,Y) 



fc±>0, k>l S^CY^ 



X 



/ ■■• / I{x,k,Oll^r Ild^^, (4) 

JCr JCr 



where k = k^ + k^. For the exponent of r we used the bilinearity of a and the obvious 
fact a{S.,Y+) = 0. 

We rephrase Lemma 5.1 of [12] to make it compatible with the present notation. 
Let A be a finite set of indices and let 

9 = 9{ii)ieA 

be a function that is analytic for all 7^ 0. Assume that foii<k 
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as C,k oo, uniformly when all the C,e with i ^ i, k are bounded and bounded away 
from zero. For B G A denote hj qb the function obtained from g by setting all C,i 
with i ^ B equal to 1. (In particular g = gA-) Define 



i<3 



Then when p, q ^ 0, 

^'-R i£A BcA ^ ieB 

where r is so small that all the zeros of the denominators lie outside Cr- When B is 
empty the integral on the right side is interpreted as . . . , 1). 

We shall apply this to the integral on the right side of (jl]). In order to do this we 
change the set of indices in that integral from [— /c_, —1] U [1, k+] to S_ U 5"+ in the 
obvious way. In the application of the lemma we take A to be the set of negative 
indices, and g the integral over the remaining variables C,i with i > (i.e., i G 5*+): 



9i0 



j>o 



Then the integral on the right side of (jl]) equals the left side of ([5]). 

Let us see that g{^) satisfies the required conditions. Since R is arbitrarily large, 
g{^) is analytic for {i < 0). When i < k < and we set = p/{l — q^k) the 
terms in the integral defining g that involve are 



As ^jt — 7- oo the product of the /s as well as the last factor are clearly bounded as 
oo. The product of the remaining two factors is 0(^^'~**'') — )■ since Sj < s^- 

So the hypothesis of the lemma holds. In its applicatior|^ we replace |v4\i?| = 
k_ — \B\ oi the with i < by 1. Each replacement has the effect of reducing the 
number of variables in g{!^) by one and dividing by so in the integral over the 
5-variables we must divide by q^+^^-~\^\\ Therefore to get the coefficients of the 



^In the resulting integrals over Cr we must have rR large, since to maintain the analyticity of 
g{£) for small (with i < 0) we must have all the (with j > 0) even larger. 
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resulting integrals over large and small contours we must multiply the coefficient in 

by q-k+{k^-\B\)_ rjj^g ^gg^i^ jg 

+ ^<r{B, S-)-k- \B\ P 



JB\{\B\+l)/2 



gk-{k--l)/2 + k+ {k--\B\) + \B\. 



Here we used the identity 

a{U,V) + a{V,U) = \U\\V\ + \UnV\ (6) 

to obtain a{B, S-) = —a{S-, B) + (A;_ + 1) \B\. Thus the integral on the right side 
of (jl]) equals 



„|B|(|B|+l)/2 
fc- + [B[^-a(5_,B) P 



^fc_(fc_-l)/2+fc+ (k_-\B\)+\B\. 
BCS- ^ 



X 



f [ Iix,\B\ + K,0 U^m^r Ud^. (7) 

When \B\ + = Q the integral is interpreted as 1. 

If we multiply the coefficient in ([7]) by the coefficient of the integral in (jlj) and use 

k{k - l)/2 = A;_(A;_ - l)/2 + - l)/2 + k.k+ 

and the bilinearity of a we get 

/ -j^\m+fc_ + |B|^m(m-l)/2-fem+cr(S_,y_\B)+(7(5+,F) 



^ ^fc+{fc+-l)/2+(fc+-l) |B| {|B|+l)/2 



A; - 1 
m — 1 



Now instead of first summing over all B G S- and then over all S- C we 
reverse the order by first fixing B and then summing over all S- satisfying B G S- G 
Y_. We change notation, so that in the end we get integrands like those in (jlj) with 
Si G S- U Sj^. To do this we interchange B and S"., so the old k_ becomes \B\ and 
the new k_ is defined to be the new l^-l- We label our indices so that Si G when 
i < and Si G S+ when i > 0. Then our formula becomes 

Py(X™(t) <X)=J2 Yl / / ^X, k, ^ H^i'' d^i (8) 

|S±|=fe± 
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where k = + and 

Cm S S = (-1)™+''- T-^i^-^)/^-k+m+a(S+,Y) ^k+{k+-l) /2+k- (k+-l) ^k- (fe-+l)/2 



\B\ + fc+ - 1 
m — 1 



X ^ (^_i^\B\^a{B,Y-\S-)-m\B\ 
Y^DB^S- 

- ^™(m-l)/2-fc+m+o-(5+,y)+A;_{A;_+l)/2 ^fc(fe-l)/2 

151 + A;. - 1" 



^ m - 1 ■ 

When /c = the integral in (|H]) is interpreted as 1. After our change of notation the 
restriction > 1 in the sum in (jl]) becomes \B\ + A;+ > 1. 

To evaluate the sum in (M) we us^ 



ucx 

|i7|=n 



{U,X) 



T 



n{n+l)/2 



\X\ 

n 



(10) 



We write the power of r in the sum in ([9]) as 
and then the sum in (191) as 



k. cr(S- ,Y.\S-)-mk. 



,^yB\S_| ^<T(B\5_,Y_\S_)-m|B\S_| 



+ A; - 1 
m — 1 



The restriction \B\ + fc+ > 1 is the same as l-SyS*-! + A; > 1. If we apply (fTOj) with 
U = B\S^ and X = Y_\S^, then the last sum may be written 



|y-\s- 



ra=0 



mn 


'\y-\s-\' 




n + k — 1 




n 


T 


m — 1 



when A; > 0. When A; = the sum begins with n = 1. (Notice that all terms of the 
sum vanish unless m < \Y_\S- \ + k.) The sum from n = is equal tc0 



(^_iyY.\S-\^\Y.\S-\{\Y.\S^\-2m+l)/2 



k-1 

m - \Y_\S. 



^The left side equals J2i<ui 



<---<u„<\X 



^tiiH hM„^ which is the coefRcient of z" in the expansion 

of nL='i(l + zt'') which, by the r-binomial theorem [9l p.26], equals EL='o r"("+i)/2z" 1^1 

^If we divide by the n = term of the series the result is the r-hypergeometric function 
24>i{a,b;c;T,c/ab) with a = r^'^-^'^-l, b = t'^ , c = 7-'=-™+i. An application of Heine's g-Gauss 
identity [5, p. 14, eq. 1.52] with q = t gives the stated equality. We "discovered" the equality by a 
Maple computation. Doron Zeilberger [T5] showed us that it follows from the g-Gauss identity. 
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When k = the condition m < |y_\S'_| + is the same as m < lYiyS*-!, and then 
the sum from n = equals zero. Therefore the sum from n = 1 equals minus the 



n = term, which is 



Using 



-1 
m — 1 



(-irr 



m — m(m— 1)/2 



+ l)/2 = a{Y^\S-, 
we see that when A; > the sum in (jO]) equals 



k - 1 

m- |yi\5_| - 1 



and so 

^ _ / i\m+|y-\5_| ^m(m-l)/2-fc+m+cr(5+,y)+cr(y_,Y_\5_)-mlY_l+A:_(A:_+l)/2 

^m, S-,S+ V / ' 



k-1 



12) 



When k = (i.e., S- = S+ = 0) the sum in (ED equals ^-m{m-i)/2 ^^ien 

m < \Y^\ and otherwise, so 0^,0,0 = 1 when m < and otherwise. Using ( fTTl) 
we find that this is exactly what we get if we substitute /c = into f|T2|) . Therefore 
^ holds for all fc± > 0. 

These are to be substituted into (IHl) to give Py(a;m(^) < 2;). To obtain the formula 
for Py(xm(i) = a;) we replace /(x, /c, ^) in ([8]) by A/(x, A;, ^), obtaining 

Py(x„(t) = X) = ^ 5^ C^,5_,5+ / / AJ(X, fc, n^'^ (13) 

Since the integral in ([8]) is interpreted as 1 for each x when A; = 0, the integral here is 
interpreted as 0. Recall that in both formulas Si G S_ when i < and Si G 5+ when 
i > 0. 



III. Formula for Py(r7^(i) = 1) 

To obtain Py(?7^(t) = 1), the probability that site x is occupied at time t given 
the initial configuration Y, we replace the coefficients in ( fT3l) by the sum of (JT2i) over 
all m. The sum of those terms involving m is 



^^(^_X)™ ^m(m-l)/2-(fc+ + |Y'_|)m 



A;- 1 
m-\Y_\S-\ - 1 
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and after making the substitution m — )■ m + iF-yS-l this becomes 

fl\\Y-\S-\^\YAS-\i\Y-\S-\~l)/2~{k+ + \Y-\)\Y-\S-\ ST^ f-^^m ^m(m-l)/2-mk ^-1 
^ ' ' m - 1 

The factor outside the sum may be written 

^-a(y-\5_,y_\5_)-feiy_\5_i^ 

while another apphcation of the r-binomial theorem shows that when A; > 1 the sum 
equals 

fc-i 

If we use 

a(r_, - (T(r_\S_, Y-\SJ) = -a{Y^\S^, S_) + (14) 

cr{S+, Y) = a{S+, Y+\S+) + k+{k+ + l)/2 + k+\Y^\S-\ + k+k., (15) 

which are obtained using the bilinearity of a and (I6j), and multiply by the factors in 
f lT2|) not involving m, we obtain for the sum over m 

k-l 

cs s+ = (-1)'' r'^('5'+.^+\-5+)-'^(^-\5-.s'-) gk(k-i)/2 Y^^^ _ ^j-^ 

i=i 

k-l 

= _^a{S+,Y+\S+)~a(Y-\S-,S-) Y^^^j _ ^jy ^^q-^ 

Thus, with these coefficients we have 



(17) 



fc±>0 S±CY± 



where the integral is interpreted as when /c = 0. 

IV. ¥{r]x{t) = 1) for Bernoulli Initial Condition 

We now derive a formula analogous to f lT7|) for Bernoulli initial condition on Z, 
with density p_ on (oo, 0] and density p+ on [1, oo). It will be a sum over k± only. 
We shall begin with Z replaced by (— M, A^], apply (ITTj) . and then let M, — > oo. 
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The probability for an initial configuration F = y_ U y+ is 

pL^-|(i-p-f^-i^-ip;'-|(i-p+)^-i^-i (18) 

To obtain the coefficient of the multiple integral corresponding to the indices k± we 
have to multiply the above by cs_,s+ as given by (fT6|) . then sum over all Y± D S±, 
and finally multiply by 11^""*' ^^^1 sum over all S± with |5'±| = k±. Because of the 
structure of cs_,s+ we need only compute the two limits 

E E p?'^'(i-p+)"^-'^^'r'^^"^'^^^"^^nc"' (19) 

S+C[1,]V] y+DS+ i>0 

is+i=fe+ y+c[i,]v] 
where 5'+ = {sj, . . . , Sk+} with si < ■ ■ ■ < Sk+, and 

JToo E P^'(i-P-)^-'^-'--^^^-\"-"-^n^r% (20) 

S_C(-M,0] y_DS_ j<0 
lS_| = fe_ i+C(-A/,0] 

where S- = . . . , with < ■ ■ ■ < s^i, multiply them together, and then 

multiply by 11^=1 (P'' ~ ^'')- fflSj) the are large and in fl20l) the are small. 

The limit f|T9l) with the factor r'^('5+'^+) instead of r'^('^+'"*^+\'S'+) was computed in 
Using that result we find that f|T9l) is equal to 

ne.---6,-i/p+(i-^''+-^+^)' ^^^^ 

We required that these be sufficiently large to take the N oo limit. 
For the limit fl2UI) we set 

Si = -s.i + 1, S+ = -5_ + 1, Y+ = -Y_ + 1, 

(observe that si < ■ ■ ■ < Sk_ ) and use 

a(yi\5_, 5_) = a(5+, yv\5+). 

If we set C,i = ^Zi, {i = 1, . . . , /i:_) then fl20|) becomes 

ji?^ E E p^-i(i-p-)*^-i^-ir--(^-^-\^-)n^~r'^'- 

S+C[1,...,A/] *>o 
S+|=fc_ Y^C[l,...,M] 



^The computation began with with formula (5) of the cited paper and the hmit is the displayed 
expression preceding identity (8). 
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The limit formula f|2T|) with an obvious modification shows that ( |20]) equals 

We require that these be sufficiently small to take the M ^ oo limit. 
Denote by ^p±{k±, p^-,^) the two limits, 

k+ 

P+ 



ip+{k+,p+,i) = Yl 



f-{k-,p-,o = n 



\^^■■■^k^-l + P+a-r'+-^+'y 

P_Cl 



A^i(e-^---e-fc_)-i-i + p_(i-r-^' 

If we recall fll6p then we see that we have shown, formally, that that for Bernoulli on 
Z with densities p± 

k-l 

nvtix) = 1) = - ^ n(^' - 

k±>0 j=l 



(22) 



The reason for the qualification "formally" is that we must justify the interchange 
of the limits M, iV — )• oo with the sum over k±. To do this it is enough that we obtain 
a convergent series when the integrands in f|T6l) with their coefficients are replaced by 
their upper bounds over M and A^. Now not only was the limit (1191) determined in 
[H], but an expression was found for finite N, namely 



n«- E n 

i=l ti>0 i=l 



l-p+(l-r^-'+^^^'' 

ii-- ■ ik+ 



This gives the N ^ oo limit if i? = \C,i\ is sufficiently large (depending on r) and also 
a uniform bound 

^kl j^-kl/2+o{k+) ^ ^23) 



for some A and arbitrarily large R. Similarly for (1201) we get a uniform bound 

j^kl ^fci/2-0(fc_) (^24) 

for some A and arbitrarily small r. 
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As for the factor AJ(x, k, ^) in the integrand, we have the easy bounds 

n \f{^^,Q\<A'lR-'l/''-<'^'^\ (25) 

0<i<j 

n \fi^um<A'-r'-/'-''''-\ (26) 

i<j<0 

ni/te,o)i<^'-'+^-'-'+- (27) 

i<0 
j>0 

(For the last we used rR S> 1; see footnote [2l) The rest of the integrand is at most 
y^k+R+k.r-^ j^o{k+) ^-o(k.) _ If we use /c_fc+ < k'i/2 + /c^/2 and combine the above 
bounds we get the bound 

jt^kl+kl+k+R+k-r-'^ j^-k%+0{k+) ^fci/2-0(fc_) 

If we take R sufficiently large and r sufficiently small (depending on this A) we get a 
convergent sum of integrals. So (!22i) is justified. 



Remark. For SSEP the only nonzero terms in fl22l) are those with A;_ = 0, /c+ = 1 
and = 0, A;_ = 1, and we get a sum of two single integrals with p± appearing only 
as coefficients: 

¥iv,{x) = l) = p+ + P- /, (28) 

Symmetrization 

In [H], where p_ = 0, we symmetrized the integrand in fl22l) . found that it was 
a determinant, and this led to a representation for W[xm{t) < x) as an integral of 
Fredholm determinants. Although it does not lead to as nice an expression, we can 
do an analogous symmetrization here. 

Write the integrand in (12^ as 

(i-n^r')n/(^-^^) (29) 

Cx s{B,')t 

X n /(^-^^■) Ht^^-^^-""-'^) 

z<i<o 2<o 

X n /(^-^^■) nT^^+^^+"'+'^^- ^^^^ 

0<i<i i>0 
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The factor fl29|) is symmetric separately in the C,i with i < and the with i > 0. 

The factor (I3T1) is a function of the .^j with i > 0. Its symmetrization, given by 
formula (9) of [H], is 

(All indices positive.) By the identity (3) of [I3] this equals 



(~^)^^ „-fc+(fc+-l)/2 1 . ^ ^ TT P+ - P) e^t 



where 



det(ir+te, 0))i<.,<fc+, 



For the factors 

i<i<0 

we write as before = ^zl, and this becomes 

TT ^ 5' ~ ^-^ TT P- 

iio p^^i^ + 9 - f=i ■ ■ ■ 4, - 1 + P-(i - r^-'-') 

k— ~ 

= TT ~ TT P- 

oi4 ^ + - i=i ■ ■ ■ 4+ - 1 + p-(i - ^^-^-^) 

The symmetrization of this is 



J_ /c_(fc_-i)/2 TT ^ TT ^ P^i^ .^^-^^^ positive) 



= -l^/-(^-^)/^ n A/^ . n ,-i K (all indices negative). 

Hence the symmetrization of (!30|) is (all indices negative) 
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kj. p lY^^'' 11 - 1 + P-(i - r-^) -^Z 1 - 



^ U + a^e. - 6 J 11 -I + PJI- r-i) 



VP + - 6/ V C - 1 + P-(l - r 



where 



e~i-l + p_(l-r-i) p + ga'-f 
So our formula may be written 

fc_+fc+-i 

'.■_,'=+ >o ■ +■ j=l 
fe_+fe_l_>l 

X L L(i-nc^)n/te'0) 

JCj^ JCr j i<0 

• n^^* ■ n^^- (32) 

j<0 i>0 

Remark 1. Suppose the change of density occurs at y rather than zero. Then 
to use the computations of (|T9l) and (!20l) we would make the substitutions Y± = 
Y± + y, S± = S± + y, Si = Si + y and apply those computations to Y±, etc. The 
result in the end is that each in the formula for A/(a;, k, ^) becomes ^f^^. In case 
P- = p+ it makes no difference what y is so the formulas we got were independent of 
X, as they should be. 

Remark 2. To obtain formulas for the correlation function F{r]t{x) = 1, rio{0) = 1) 
we write it as 

¥{7]t{x) = 1) - ¥{7]t{x) = 1, r/o(0) = 0). (33) 



The first probability we know. For the second we modify (U81) by multiplying by 
1 — p_, the probability that site is initially unoccupied, and take only those YL that 
are contained in (— cxd,— 1]. To use the preceding computation we use the one-one 
correspondence between subsets S- C (— oo, — 1] and subsets S'_ C (— cxd,0] given 
by S'_ = + 1. Then a{Y-\S-, S-) = a{Y!^\S'_^, S'_) and, with obvious notation, 
Yl ^j"'*' = n ^ n 4 • follows that for the second probability in (l33l) we multiply 
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the integrands in ( 122|) and fl32l) by (1 — p_) ni<o'^«- Therefore for (l33l) itself we 
multiply the integrands by 

i-(i-p-)n^- 

i<0 

V. Formula for the Total Flux 



Suppose maxF_ < < miny4.. When t = there are |yi| particles < 0, so Qt, 
the total flux to the left across at time t, is the number of particles < at time 
t minus Thus Py(Qt > m) = Py(x^+|y_|(t) < 0). Therefore in (ED and 

we replace m hj m + \Y_\ and set x = 0. If we use (fT^ and ( [T5|) and replace m by 
m + (so the new m may be negative) then f ll2p becomes 

_-|^-jm+A.- ^o-{S+,y+\S+)-(T(y_\S-,5_)+m(m-l)/2+fc+(A;++l)/2-mA:+ 



X g 



k{k-l)/2 



Cm+\Y-\,S-,S+ - 

' k-1 ■ 

m + — 1 

Using this, comparing with f|T6|) . and following the argument of the last section we 
obtain the formulas 

k-1 



¥{Qt > m) = (-1)"+^- ^m(m-l)/2+A:+(fc++l)/2-mfc+ gfc(fc-l)/2 

fc±>0 



J Cm Cr 



(34) 



fc±>0 



A;_! A;, ! 



A;- 1 
m + A;_ — 1 



X 



/ fc+ / ri'^*^^^' '^j) ■ '^i))-fc-<i,i<-i • det(K+(^i, Cj))i<i,j<fc+ ■ n^^i 



j>0 



where in the expressions for K±{^i, ^j) we set x = 0. 

For P(Qt = m) we subtract from these what we get by replacing m by m + 1. To 
get a formula for (e'^'^*), the expected value of c*^*^*, we then multiply by e^"^ and 
sum over all m. 



If ¥{Qt >m)=pm then (e^'^'^) = Z]m=-oo(2m ~ Pm+i) e^"". The series converges 
for all A and represents an entire function of A|j Therefore if we find a formula that 



^For positive m, if Qt > m then the particle initiahy the mth to the right of must have moved 
at least m steps to the left at time t. The probability of this is less than the probability that 
a free particle would so move, and this probability is o(e-" '°s'"+<^(™)) as m -> +oo. Similarly 
]t> m) = l~ 0(e-l'"l i°g|™l+o(lm|)) as m -oo. 
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holds when 3? A > 0, and the formula represents an entire function of A, it will hold 
in general. When 3? A > we can write the series as the difference of two convergent 
series in the obvious way. Changing the summation variable in the second series and 
combining gives (1 — e~^) J2m=-ooPm^^"^- This is what we do below. 

The factor that involves m in the first formula is 



^ ]_^™' ^»7i(r?i— 1)/2— m k+ 



k-1 

m + — 1 



When we make the replacement, subtract, multiply by e , and sum over m we get 



k+ 

m=—k^ +1 



m ^m(m— 1)/2— m A;+ m 



-l)'"r 



m + fc- — 1 



and another application of the r-binomial theorem shows that this equals 

fc-i 



\k+ 



1 



g-AN gAfc+ ^-fc+(fc+ + l)/2 



nil 

i=i 



k~l 

-fc+(fc++l)/2 J-j-j^-^ _ 



Therefore 



fc±>0 



k ^k(k-l)/2) gAfc+ 



i=o 



X 



JC'1+ JCr' 



To justify the interchange of the sums over m and k±, and to show that the 
sum represents an entire function of A, we use estimates analogous to those used in 
Sec. IV. For the integral dMD we have the bound A''^+''+^+''-''~' R-kl+o{k+) ^kl/2-o(k.) 
and for its coefficient another bound for all m, since \m\ < k when the coefficient 



is nonzero. And |e 



Am I 



< e 



for such m. If i? is large enough and r small enough 



the sum (over k± and m) of the bounds is finite and uniformly bounded in bounded 
A-sets. This justifies the computation. 

In the case of SSEP (r = 1) our formula simplifies since 



and we can write it as 
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fc±>0 

^ /(o, k, ■ f-ik-, 1, f+{k+, 1, n • 

In particular, this is a function of the two quantities p_ (1 — e~^) and {e^ — 1). 

This fact was derived in |3j in a different way. The authors went further to show 
that it is a function of a particular combination of these quantities. This allowed 
them to reduce to the case when p_ = and to obtain, for general p±, a series for 
the generating function involving integrals only over small contours. 
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